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Abstract 

Studies on the isospin of fragments resulting from the disassembly of highly 
excited large thermal-like nuclear emitting sources, formed in the Au+ 197 Au 
reaction at 35 MeV/nucleon beam energy, are presented. Two different decay 
systems (the quasiprojectile formed in midperipheral reactions and the unique 
source coming from the incomplete fusion of projectile and target in the most 
central collisions) were considered; these emitting sources have the same ini- 
tial N/Z ratio and excitation energy (E* ~ 5-6 MeV/nucleon), but different 
size. Their charge yields and isotopic content of the fragments show different 
distributions. It is observed that the neutron content of intermediate mass 
fragments increases with the size of the source. These evidences are consistent 
with chemical equilibrium reached in the systems. This fact is confirmed by 
the analysis with the statistical multifragmentation model. 
PACS numbers: 25.70Pq, 25.70-z, 24.60-k 



The study of heavy-ion collisions at intermediate energies (10<E<100 MeV/nucleon) is 
a useful tool to investigate the mechanisms of fragment production in highly excited nuclear 
systems. In this energy regime multifragmentation appears as one of the main de-excitation 
channel |IJ, which can be considered 0,f| as a manifestation of a liquid-gas type phase 
transition in finite nuclear systems. This gives an access to the nuclear equation of state. It 
has been shown that statistical models 0|| happen to be very effective in the reproduction 
of main characteristics of the fragment production such as charge distributions and fragment 
correlations. 

Nowadays, with the advent of radioactive beam facilities, the influence of the isospin de- 
gree of freedom is strongly addressed and experimental information on the isotopic content of 
emitted fragments represents a meaningful starting point in order to get a deep understanding 
of either the de-excitation mechanisms or the nuclear matter properties. 

Theoretical calculations predict that, in the fragmentation of asymmetric nuclear matter, 
the isospin composition of the "liquid" phase (usually associated with intermediate mass 
fragments (IMF) and heavy residues) and the "gas" phase (light charged particles (LCP) 
and nucleons) depends on many factors. In particular, calculations by Miiller and Serot IJ 
showed that, for very neutron-rich systems, there may exist a distribution of the excited 
nuclear matter into a neutron-rich gas and a more symmetric liquid. However, Lee and 
Mekjian |7| have pointed out that the Coulomb and surface effects, which are important 
for finite systems, may moderate the neutron enrichment of the gas phase by producing more 
free protons. The Statistical Multifragmentation Model (SMM) calculations || have shown 
that the neutron content of the IMF can increase in the region of the phase transition, in 
agreement with the experimental observations M. 

In this situation only new experimental data could solve the problem of isospin composi- 
tion of the gas and liquid phases. Indeed, a variety of experiments can be found in literature. 
In Refs. |T0|JTT| the difference of the mean N/Z ratio of the LCP and of the IMF has been 



interpreted as the separation of gas and liquid phases. The authors of Ref. [0 considered 
the presence of neutron rich LPC and light IMF as an evidence of the neutron enrichment of 
the gas phase. The average N/Z ratio of fragments emitted from excited nuclear systems is 



seen to vary with the excitation energy and the N/Z ratio of the system. In particular in the 



work of Ramakrishnan et al. [fL3|1 , where for different reaction the beam energy and the mass 
of the system were kept constant, it has been shown that the IMF's isospin present a linear 
dependence on the N/Z ratio of the system. Moreover, it has been observed that the neutron 
content of the emitted IMF depends on the excitation energy of the (fixed size) sources ||, 
increasing as the excitation energy increase. 

In this paper we present experimental data extracted from two emitting sources with 
approximately the same N/Z ratio and excitation energy, but with different size. Since existing 
experimental data cover the study of the N/Z ratio of fragments as a function of the excitation 
energy and N/Z ratio of the emitting source, this analysis represent a complementary point 
to fill the experimental picture on the effects induced by the isospin on the decay process. 

It is also important to note that these new data provide information about isospin of 
fragments produced in decay of the largest nuclear systems under investigation up to now 
(A>300). This allows for more reliable extrapolations for the case of nuclear matter, as well 
as for the astrophysical applications in supernovae explosions and neutron stars. 

The reaction Au+Au at 35 MeV/nucleon was studied in experiments performed at the Na- 
tional Superconducting K1200 Cyclotron Laboratory of the Michigan State University. Light 
charged particles and fragments with charge up to Z=20 were detected at 23° < 9i ab <160° 



by the phoswich detectors of the MSU Miniball hodoscope [T3j. The MULT ICS array [15] 
covered the angular range 3° < Q\ a b <23° and allowed for a charge discrimination up to 
Z=83 and for good mass discrimination for Z=l-6 isotopes. The geometric acceptance of the 
combined array was greater than 87% of An. 

A brief summary of the results published in so far will follow. The selection of the impact 



parameter b is based on the number N c of charged particles detected [ |16 

1/2 



b/b max = P(N' c )dN')j 



where P(N C ) is the charged particle probability distribution and tt ■ b 2 max is the measured 
reaction cross section for N c >3. 

Through the analysis of energy and angular distributions of the experimental data, it was 
possible to identify events originating either from disassembly of a unique source formed in 



central collisions, or from decay of the quasiprojectile (QP) in peripheral and mid-peripheral 
collisions ||17|| . The excitation energy of the fragment sources was calculated both via the 
calorimetric method and through comparison with model calculations Hl7||l8| . Extra kinetic 
energy, such as radial flow or rotational motion, does not contribute appreciably to the 
excitation energy. This results from the analysis of charge, angular and kinetic energy distri- 
butions of the emitted fragments, as well as from the study of event by event charge partition 
17|- |i9pT| . Finally, nuclear temperatures were determined using the technique of the double 



ratio of isotope yields |20|j . 

The analysis determined that the thermal characteristics (Excitation energy E* and 
isotope temperature T iso ) of the QP, formed at 0.6< b <0.7 (E*=5.5±0.6 MeV/nucleon, 
T iso =4.3±0.2 MeV), and of the unique source (E*=5.5±0.5 MeV/nucleon, T iso =4.3±0.4 
MeV) are remarkably similar [|l7j. 

It is therefore important to analyse the whole picture of fragment production in these 
reactions. The aim of this paper is to provide new results by comparing IMF production 
between emitting sources of similar excitation energies (E* ~5-6 MeV/nucleon) and N/Z 
ratio (118/79=1.49), but of different size (Z=79 for the midperipheral QP, Z=126 for the 
central unique source (CUS) [fLgfl). In particular we will focus our attention to the isospin of 
emitted fragments. 

Charge yields of fragments obtained for the Au+Au reactions at 35 MeV/nucleon are 
presented in Fig.l. The heaviest source (CUS) decays emitting fragments lighter than the 
QP, and in a number larger than what expected by a simple scaling factor. In other words, 
the CUS undergoes a stronger disintegration at the de-excitation stage. 

The different partition in the two systems affects the number of neutrons in the IMF. 
Fig. 2 shows the isotope production yields for different Z values. The IMF emitted from the 
CUS are more neutron rich than those from the QP. This effect is enhanced when plotting 
the relative yield ratio vs the neutron excess (see Fig.3a): the ratio between the CUS and QP 
yields increases with (N-Z). In Fig. 3b we present the ratio of relative yields of neutron-rich to 
neutron-poor isotopes at fixed Z values, for both CUS and QP. It appears that the CUS emits 



preferentially the more neutron rich fragments. The average <N>/Z value of each atomic 
specie versus its charge Z is presented in Fig.4. 

One can clearly see from all the figures that the disassembly of the CUS system into 
a relatively larger number of fragments leads to production of more neutron rich IMF and 
LCP. In the previous study it was experimentally demonstrated that fixing the size of the 
source, but varying its excitation energy, the increase in the multiplicity of emitted fragments 
is accompanied to more neutron rich IMF ||. In the present case we observe that this effect 
is even more pronounced. 

These results are consistent with the statistical picture of disintegration of finite nuclear 
systems. For example, the statistical multifragmentation model (SMM), see e.g. ||, is based 
upon the assumption of statistical equilibrium at a low-density freeze-out stage. Different 
break-up partitions are sampled, according to their statistical weights, in the phase space. 
After break-up of the nuclear source the fragments propagate independently in their mutual 
Coulomb field and undergo secondary decays. The de-excitation of the hot primary frag- 
ments proceeds via evaporation, fission or Fermi-break-up p2 |. SMM is very successful in 
reproducing experimental data concerning both peripheral and central nucleus-nucleus col- 
lisions |23^p5|. In particular a very good agreement between experimental data and SMM 



calculations was found in the study of the Au+Au 35 MeV/nucleon reaction [17-19|. Here we 
use a new version of SMM based on the generation of a Markov chain of partitions [H] . This 
version keeps its full reliability concerning the charge distribution predictions and allows for 
taking into account all effects influencing the isotope content of the produced fragments. 

In this paper, by basing on this statistical model we aim to single out the main physical 
effects responsible for the observed trend. Even though some assumptions of the model, 
such as a fixed freeze-out volume or non-overlapping fragments, are approximations of the 
real conditions, and some parameters (flow energy, fragment's level density) are not clearly 
defined, these uncertainties influence the fragment charge distributions mainly and can be 
accounted for by finding the source distribution with the well elaborated technique f2l |. When 
the charge distributions are fixed, the main features of the studied isotope distributions are 
driven by the chemical equilibrium effects and the binding energy of the fragments. 



The events generated by SMM were filtered to take into account the experimental effi- 
ciency. In the Figures 1, 2 and 4 the comparisons between experimental data and SMM predic- 
tions is shown. The following set of parameters gives a good agreement between experimental 
data and calculations (see also 18 ]2jJ ) : an excitation energy i?*=5.5±0.6 MeV/nucleon; 
A s =197, Z s =79, p s = |p , for the QP; A s =315, Z s =126, p s = \p Q for the CUS (p =0.15 fm~ 3 
is the normal nuclear density). The smaller density in the central collisions is consistent with 
an additional expansion caused by the flow development. The microcanonical temperatures 
obtained in both cases are also quite similar: T m i cr « 5.4 MeV for the central case and ~ 
5.5 MeV for the peripheral one. 

Dynamical calculations P3,2E,Z7| predict that N/Z ratios of the emitting sources must 
have values close to that of the initial system. In particular an analysis dedicated to the 
study of 50 MeV/nucleon central collisions of 112 Sn+ 112,124 Sn f2?J has shown that the best 
agreement between calculations and experimental data claims for N/Z values of the source 
close to that of the starting system (differences are lower than 3%). Also the dynamical study 
of 129 Xe+ 197 Au central and peripheral collisions at 50 MeV/nucleon has shown that the N/Z 
ratio of the thermal sources does not change essentially (differences are within 2%) from the 
initial one [p6fl . Results presented in Ref. |23) and concerning the study of peripheral Au+Cu 
reaction come to the same conclusions. Generally, since the considered emitting sources are 
large in size, fluctuations in their N/Z ratio value should be strongly reduced. Therefore, the 
assumption of conservation of the N/Z ratio after the dynamical stage, seems quite realistic. 
We note, that our following interpretation will be valid even under less strict assumptions, 
namely, when the N/Z ratio changes in the same way in the both cases, or it becomes smaller 
for the CUS. 

The reason of why the nuclear sources at practically the same temperature and excitation 
energy produce so different fragment charge distributions (Fig.l), can be found in the larger 
Coulomb energy in the CUS case, and in the different dynamics of formation of the CUS 
leading to a flow and, as consequence, to smaller freeze-out densities. 

One can see also from the Figures that the calculations reproduce the trend of increasing 
neutron content of produced fragments with disintegration of the nuclei into smaller pieces. 



This trend can be explained as follows: if the chemical equilibrium is established the big 
fragments have larger N/Z ratio than small ones. After disintegration of the big fragments, 
independently if caused by increasing excitation energy (or temperature) of a thermal source 
as in H or by a size (Coulomb) effect, the neutrons of the big fragments are accumulated 
mainly in the small secondary fragments and not in free neutron gas. The SMM calculations 
predict that the mean number of primary free neutron in the freeze-out volume increases 
only from 0.98 for the peripheral source to 1.79 for the central source. This gives rise to the 
N/Z ratio of the IMF and LCP, which is preserved after the secondary de-excitation. 

Even if the model well explains the observed effects, the experimental rise in neutron 
content is slightly more pronounced. One can speculate, that the initial CUS could have 
slightly larger N/Z ratio than the peripheral source. Our calculations show that a slight 
decrease of this ratio to N/Z=1.43 without changing other SMM parameters, for the QP, 
would be sufficient to explain this disagreement. However, as we pointed above the difference 
(if any) of this ratio of the sources should not be significant. For this reason, we believe that 
the predicted redistribution of neutrons from heavy fragments to light ones with disintegration 
of a nuclear system is the natural mechanism for explanation of the data. 

In summary in the study of the 197 Au+ 197 Au 35 MeV/nucleon reaction it was possible to 
well identify two different emitting sources with the same N/Z ratio, excitation energy and 
temperature (E* ~5-6 MeV/nucleon; T iso ~ 4.3 MeV from the experimental measurements, 
T m i cr ~ 5.4-5.5 MeV from the SMM predictions), but different size: the QP formed in 
midperipheral reactions and the CUS coming from the incomplete fusion of projectile and 
target in the most central collisions. Charge yields and isotopic content of the fragments show 
different distributions in the two cases. The CUS decays emitting lighter fragments than the 
QP and this fact reflects on the multiplicity of fragments from CUS, which is higher than 
what expected scaling for the size factor the QP multiplicity. Thus, even if some parameters 
are fixed (E*, T, N/Z) the partition of the system depends on the size and the density which 
determine the Coulomb energy in the freeze-out. Moreover, the neutron content of light 
charged particles and intermediate mass fragments increases with the size of the source. 



Primary fragments produced in the freeze-out are hot and decay afterwards by emitting 
mainly neutrons and light charged particles. Therefore, it is natural to connect the behaviour 
of the experimental isotope yield as a function of source size with the similar evolution of 
the N/Z ratio of the corresponding hot fragments. This conclusion is here supported by the 
SMM calculations which reproduce the fragment production in a reasonably good way. The 
data here presented, obtained both in central and peripheral |J collisions, indicate that the 
neutron content of the fragments, produced by the decay of thermal-like systems, increases 
with the multiplicity of emitted IMF. The data are consistent with hypothesis of thermal and 
chemical equilibration in finite nuclear systems ||, which leads to production of IMF with 
large neutron content. This justifies the method of extracting temperatures through isotope 
thermometers [17] and the thermodynamical description of these reactions. 
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new version of SMM. The authors would like also thanks J.D.Dinius, S.Gaff, C.K.Gelbke, 
T.Glasmacher, M.J.Huang, W.J. Lynch, C.P.Montoya, M.B.Tsang and H.Xi, for their collab- 
oration during measurements and data analysis. One of the author (A.S.B.) appreciate warm 
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Figure captions 



Fig.l: Charge yields of fragments obtained for the Au+Au reactions at 35 MeV/nucleon 
fl7|J18ll . Circles represent the experimental data: the solid ones are for the QP systems, the 
open ones are for the CUS. Solid and dashed lines are SMM calculations for peripheral and 
central cases, respectively 

Fig. 2: Relative yields of different isotopes for fragments with charges from Z=l to Z=6. 
Notations are as in Fig.l. 

Fig. 3: (a) ratio between the CUS and QP yields for each isotope as a function of N-Z. 
(b) ratio of relative yields of neutron-rich to neutron-poor isotopes at fixed Z values, for both 
cases (solid point are for the QP, open points for the CUS). 

Fig. 4: Mean neutron-to-proton (< N > /Z) ratio of the produced fragments. Notations 
are as in Fig.l. 
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